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Abstract 



A theory of a lateral stress relaxation in a fluid bilayer membrane 
under a step-like pressure pulse is proposed. It is shown theoretically 
that transfer of lipid molecules into a strained region may lead to a 
substantial decrease of the membrane's free energy due to local relax- 
ation of the stress. Simultaneously, this same effect also causes ap- 
pearance of the spontaneous curvature of the membrane. Proposed 
stress relaxation mechanism may explain! recent experimental ob- 
servations of a collapse of the ionic conductance through the protein 
mechanosensitive channels piercing the E-coli's membrane!'!. The 
conductance decreases within seconds after opening of the channels 
by applied external pressure. The time necessary for a propagation 
of a stress over the largest membrane's dimension is 5 -r- 6 orders 
of magnitude shorter. On the other hand, theoretical estimates of 
the phospholipids diffusion time inside the strained area is favorably 
comparable with the experimental data!. Our theory also predicts 
that the membrane's curvature increases simultaneously with the 
stress relaxation caused by the lipids diffusion. This effect is pro- 
posed for experimental check up of the theoretical predictions. 
82.65.-i, 68.15.+e, 46.30.-i 



I. INTRODUCTION 

The functioning of the protein channels (pores) in the cell membranes reg- 
ulates flows of ions in- and out of the cell, thus influencing signal transmission 
processes in the neural networks of the biological systems^. Therefore, the gat- 
ing mechanisms of the channels are of great interest for the fundamental and 
practical matters, e.g. like enhancement of the efficiency of the general anesthe- 
sia methods!. Mechanical stresses in cell membranes may cause opening of the 
protein mechanosensitive (MS) channels!, which hence play a role of mechanore- 
ceptors in the cellular organisms. One of the biological objects convenient for 
experimental study is the large conductance MS channel (MscL) in the inner 
membrane of the Escherichia coli (E-coli) bacteria. It is possible to measure a 
single channel conductance, an find its dependence on the internal lateral ten- 
sion in the membrane. The tension, which causes opening of the channel, could 
be evaluated using video microscopy measurements of the membrane curvature 
formed under the applied external pressure!. 

The time- dependence of the membrane's conductance observed by C.C. Hase 
et a/.! reveals gradual collapse of the ionic currents through the MscL's within 
seconds after their opening. According to existing hypothesis!, this may happen 
due to mutual slide of the lipid monolayers, constituting a bilayer membrane of 
the E-coli. This slide would then encourage a relaxation of the lateral stress inside 
the membrane, which is necessary for an opening of the MscL's. The slide of the 



monolayers should be encouraged by the same pressure pulse, which causes initial 
opening of the ionic channels. 

In this paper we propose a theory of such a stress-relaxation process and 
demonstrate that the slide of the monolayers under a constant external pressure 
is indeed thermodynamically advantageous. In order to show this, we calculate 
the free energy difference between two possible equilibrium states of a membrane 
under a pressure gradient pulse. 

The first ("intermediate") equilibrium state, is achievable shortly after the 
beginning of a pressure pulse P, within the time of propagation of the mechanical 
stress along a membrane At > 10~ 6 sec (evaluated below). In this state the initial 
numbers of the phospholipid molecules in the monolayers constituting the bilayer 
membrane, iV\ and N 2 , yet, remain unchanged: Ni = N 2 = N . Here N is 
the number of the molecules in a monolayer before external pressure gradient is 
applied. We neglect a small probability^! of the flip-flop process, i.e. exchange 
of lipids between the monolayers. Then, the tensions T\, T 2 , in both monolayers 
approximately coincide, since the radius R, of the membrane is of order 1 -i-2fim, 
while its thickness Iq, is about 2nmi. 

In the second, (final) equilibrium state, achievable within diffusion time r ~ 
lsec 3> At, the extra lipid molecules from the outer patch of E-coli's material 
have enough time to be "sucked in" inside the concave-side monolayer, see Fig. 1, 
until the tension in the monolayer is relaxed at a finite curvature: T\ = 0. Then, 
the tension T 2 in the convex-side monolayer (which is stuck to the pipette and 
thus can not absorb extra molecule JUl, see Fig. 1) alone maintains a mechanical 
equilibrium under the same external pressure difference P across the membrane. 
In accord with the Laplace's law we have: 



Here r is the radius of curvature. Next, we calculate self-consistently the free en- 
ergy of the whole system: membrane+patch, in both equilibrium states. Impor- 
tant, the free energy in the relaxed, second state, proves to be lower than in the in- 
termediate state with constrained numbers of lipids, Fig. 2. Choosing the second, 
relaxed state we then evaluate the number AiV\ of the " sucked in" lipids and find, 
see Fig. 4, that it is of a macroscopic magnitude, e.g. ANi ~ 15%iVi. Based on 
these results and using experimentally determined diffusion coefficient of individ- 
ual lipid molecules in bio-membranes at room temperature! D ~ 1 4- 20/im 2 /sec, 
we estimate the characteristic diffusion time r: 



which is comparable by the order of magnitude with the experimental dataB. 
Thus, we predict a macroscopic increase of the number of lipid molecules AN± in 
the concave monolayer and, hence, propose to measure the concomitant change 
of the weight of the part of the membrane sucked in in the .pipette during the 
decrease of the ionic conductance (in the setup of the worksS!). 



T 1 + T 2 = T 2 = Pr/2; at T x = 0. 




AN±R 2 
iVi ~D 



~ lsec, 



(2) 



Besides, as it is apparent from Fig. 2, the equilibrium curvature of the mem- 
brane (indicated with the arrows) also changes significantly as the membrane 
passes from the intermediate to the final equilibrium state. Hence, it should be 
possible to check our theory also by measuring the time-dependent increase of 
the curvature, which should then accompany the decay of the ionic conductance 
observed! under a constant external pressure gradient. 

The plan of the article is as follows. In Section II we describe a free energy 
functional in our model of a bilayer membrane and derive the self-consistent 
equations for the membrane in the equilibrium state under a constant pressure 
gradient. Both versions of the equilibrium state discussed above are considered. 
In Section III we give analytical and numerical solutions of the derived equations 
and discuss in detail corresponding calculated dependencies. Finally, we discuss 
in the last Section IV possible experimental verifications and future improvements 
of the theory. 



II. THE MODEL FREE ENERGY 

We use the microscopic model of a bilayer membrane as described e.g. inSi. 
First write a single i-th monolayer (i = 1, 2) free energy per lipid molecule fi. 

C k 

fi = fsi + fu + fu = l<k-\ ^irih-h) 2 (3) 

a,i 2 

Here f S i is the external surface energy of the monolayer, which increases with 
the area per lipid molecule due to the hydrocarbon tails interaction energy 
with the solvent. The short-range repulsion between the polar molecular heads 
at the surface of the membrane is represented by the second term in Eq. (|3|) . The 
third term, though looks like a spring elastic energy (with not stretched spring 
length l s ), codes for the entropic repulsion between the hydrocarbon tails in the 
depth of each monolayer^. In what follows we omit the second term in Eq.([3]) on 
the empirical groundsel relevant for the chain molecules in phospholipid bilayers. 
There is no interface energy term included in (|3|) for the internal surfaces of the 
monolayers forming the bilayer, as the latter are not reachable for the solvent 
molecules. We also assume a vanishing interdigitation of the tails between the 
adjoint monolayers. The length of a lipid molecule Zj, i.e. the thickness of the 
i-th monolayer, is not an independent variable from the aria per molecule a^. The 
volume v per lipid molecule is conserved, i.e. the latter could be considered as 
incompressibleBI. Due to this conservation condition the length U is related to the 
mean and gaussian curvatures H and K at the interface between the monolayers 
by a well known differential geometry formulas^ 

= + f&# (4) 

h = k2 + ll 2 H + hl 2 K (5) 

Here the change of sign in front of if -term in the equation (0) relative to that in 
equation (H) is due to a simple geometrical fact that the external normal vector 



to e.g. the layer 1 is simultaneously an internal normal vector to the layer 2 at 
their mutual interface. Different Zo» parameters just reflect the incompressibility 
of the lipid molecules mentioned above: 
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Substituting Eqs. ©,© into the third term of Eq. ([3D one finds contributions 
to the free enrgy from the tails in the form: 
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(8) 



where parameters have meaning of the local spontaneous curvatures of the 
i-th layer: 
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(9) 



We readily recognize the Helfrisch's formula0 for the free energy of a mem- 
brane in Eqs. (0); ©• I n a symmetric bilayer case iV\ = N 2 (where again Ni 
and N 2 signify the numbers of lipid molecules in the 1-st and 2-nd monolayer re- 
spectively) , with equivalent conditions on the opposite surfaces of the membrane 
di = a 2 , the linear in H terms in Eqs. (|7p, (||) would cancel in the total free 
energy F: 



F = E Fi = E Nifi 



(10) 



i=l,2 



i=l,2 



Nevertheless, as shown below, a linear in H term may arise when the up vs 
down symmetry of a membrane is broken, e.g. by applied pressure gradient in 
combination with the different boundary conditions at the monolayers peripheries, 
see Fig. 1. 

To make the whole idea transparent we restrict our present derivation to the 
case of a spherically "homogeneous distributions" of molecules, i.e. considering 
ctj's as being different for the different indices i's, but position-independent within 
i-th monolayer. Below, we neglect the inhomogeneous distribution of the strain 
across the thickness of each monolayer. Also, we consider only spherical shapes 
of the deformed bilayer membrane, while introducing position independent (over 
the membrane's surface) mean and gaussian curvatures: 



H = -; K 

r 



[ID 



where r is the radius of curvature. The radius of the base of the curved membrane 
is fixed at R < r, in accord with the fixed radius of the pipette, which sucks in 
the membrane, and creates a pressure difference P between inside and outside 



surfaces of the membrane in the experimental setup&tl According to experimental 
data!'! the membrane under consideration is rather thin: 



L ~ 1Q~ S R < R < r. 



(12) 



Choose direction "up" inside the pipette, perpendicular to the plane of the 
initially flat membrane. Also, ascribe indices "1" and "2" to the parameters of 
the lower and upper monolayer of the membrane respectively. Then, the external 
surface areas S\ and Si of the lower (concave) and upper (convex) monolayers 
and the numbers of the molecules in them are related by the following equations: 



Si = 2vr(r - hf 1 



S 2 = 2vr(r + l 2 f 1 



\ 
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7iR 2 = N a , at r = oo. 



(13) 

(14) 
(15) 



Here the thicknesses k=i j2 of the i-th monolayer are defined in Eqs. (|J), (§). 
A flat, undeformed membrane has a monolayer thickness Iq = v/ao under zero 
external pressure gradient, where the area per molecule ao is determined from 
the condition of the minimum of the free energy F defined in Eq. (|T0|): 



5F 

— — = , at: ai = a 2 = a , H = K = 0, N\ 
da 



No = Nn. 



(16) 



Allowing for the experimental setupB, it is reasonable to suppose that the 
number of molecules in the upper monolayer A^, which "sticks" by its peripheral 
to the pipette walls, is fixed during the experiment: 



N* = N, 



o- 



Hence, Eq. (0) together with Eqs. flT|), (|TT 
radius r to the area per molecule a 2 : 

a 2 = a 2 (H) 



(17) 

and (H) relates curvature 

(18) 



In order to find the equilibrium areas a% and a 2 one has to solve different sets 
of self-consistency equations, depending on the conservation/nonconservation of 
the number of molecules N\ in the 1-st monolayer. 

When Ni is conserved we have: 



N 1 = N = N 2 . 
This equation, together with Eqs. (0), (fLiP , 



and 



(19) 

relates a\ to the 



curvature H similar to Eq. (18 



ai = ai(H) 



(20) 



Hence, the only independent variables left are the curvature radius r, and the 
pressure difference across e.g. 1-st monolayer Pi. The equilibrium values of these 
parameters minimize the free energy in the equilibrium state of the membrane 
exposed to the total applied pressure difference P: 

5 F Pr P 

-P = T* = ^f; Pi + P 2 = P; Pi«P 2 = £. (21) 

ddi 2 2 

Here the tensions (force per unit length of the peripheral) acting in the mono- 
layers are related to the corresponding pressure differences via the Laplace's 
law (|l|). A distribution of the pressure differences Pi, acting across each mono- 
layer, generally should be determined self-consistently by a minimization of the 
membrane's free energy F, but for a thin membrane (see fljjj)) we choose them 
approximately equal. 

In the case, when N% is not fixed and the 1-st monolayer is free to absorb 
extra lipid molecules from the patch/reservoir through its peripheral edge, we 
substitute equilibrium conditions (^) by the different ones: 

SF 

f- = Ti = ; (22) 
5F Pr 

r T -Y i (23) 

P = F + ^patch • N 1+ N. 1 + Af patch = const (24) 

where for simplicity the patch free energy Pp a t c h is assumed to be equal to 
the energy of a flat piece of the bilayer membrane under zero pressure and zero 
tension, which contains N^^^ lipid molecules in total: 

F patch = E T^patch/o; fo = fi = h at: a x = a 2 = a , H = K = 0. (25) 

i=l,2 Z 

The free energies fi are defined in (||). Finally, we solve equations (|2l| ) and (22), 



(p3|). As a result, we determine the final state with the lowest energy by a direct 
comparison of the values of the free energy F in the two different cases of the 
"patch-disconnected" and "patch-connected" membrane. 



III. SOLUTIONS 



It is convenient for calculational purposes to rewrite expressions obtained in 
the previous Section using dimensionless parameters instead of a, and H. Because 
the expressions for the free energy (|T0D, ( p5|) will be studied mainly numerically, 
we shall incorporate tensions Tj into the free energy and look for the minima of 
the latter on the proper manyfold of the independent variables, rather than solve 
directly equations (|2l|) - (|23| ) of the Euler-Lagrange type. For this purpose we 
pass from the free energies fi to the new ones /j, which are Helmholtz free energies 
per lipid molecule in the i-th monolayer under an external pressure difference Pf 



rZi(h) p v 

f i = Mz i {h),h)-j iio 2hb) dz ; (26) 

A(zi, /i) = e^i + e t |y (7zr 4 - 4z7/ 3 ) + 2/i(z7/ 2 - zf 3 ) + 1 + z? - 2z^ ; (27) 

/o = /1O0, h = 0) = e s z + e t {l + zf - 2z 1 }; (28) 

f2(z 2 ,h) = f 1 (z 2 ,-h). (29) 

Here dimensionless parameters Zj and h, as well as the other ones, are defined as 
follows: 



ads h h _ rr, h Is iv k s l 2 s 

; h = - = HI,; h = — ; e s = — ; e t = — . (30) 

i> r it I ' 



A dependence hi(z) on the right- hand- side of Eq. fl2"q) is just the inverse of either 
Eq. ( PU[ ) or Eq. QT8| ) written in the dimensionless units fl3D|). It is useful for 
the comprehension of the rest of the paper to mention here orders of magnitude 
of the main parameters. As it follows from fll2]): h ~ 10~ 3 <C 1. Next, it 
follows from the definition (^) and Eqs. (f|), (|) that t> ~ Z s aj, which leads to 
an estimate: Zi ~ 1. The entropic nature of the term in Eq. (|3|) is reflected 
in the temperature dependent coefficient k s , which for the known phospholipid 
bilayers is of the orderO: k s l 2 s ~ 1 4- 10/cT, where /c is the Boltzman's constant 
and T is the absolute temperature. The scale of the surface term / s j in Eq. ( 
is characterized by the coefficient 7, which at room temperature is of the orde 
7 ~ O.lkT/A 2 . Gathering all the estimates, and also taking into account that 
the typical value of an area per molecule isi: cij ~ 60 A 2 , we obtain the following 
list of the estimates (taken for the room temperature T ~ 3007^): 

/i~10- 3 ; e s ~10- 13 erg; e t ~ 1 - 10e s . (31) 

An additional important dimensionless parameter relates characteristic experi- 
mental value of the pressure difference P ~ 50mm Hgi'i with the "microscopical 
bending energy" e t : 

— ~ 4 • 10" 4 at: e t ~ 10~ 13 erg; v ~ 1.2 • 10 3 A 3 , (32) 

where we use an estimate Z s ~ 20A for the typical length of a phospholipid 
molecule in a bilayer membrane!. 



Now, coming to a description of the results, we first solve equation for zq: 

e s + 2t t {zf - zf} = (33) 



dz 

or in the Cardano's form: 



Zq + pz + q = 0; where: p = —q = — - (34) 

Co 



and: ii)\(zy s (±)' + m' >0 . ( 35) 



The inequality in (|35| ) proves that there is a unique real root of the cubic equation 
(|34|). In the limit p> 1 the real root of Eq. ([34]), first order accurate in is: 

z « 1 - - = 1 - p-. (36) 



A. Patch-disconnected case 



In the patch-disconnected case, or during a short enough time t: At < t <C r, 
after an application of a pressure difference P across the membrane's plane, the 
numbers of the lipid molecules in the both monolayers of the bilayer membrane 
remain equal: N± = N 2 = N . Hence, the areas per molecule 0^=1,2 also remain 
equal to each other, while simultaneously undergoing an increment with respect 
to the initial undeformed value ao = vzq/I s . Neglecting small corrections oc k/r ~ 
1(T 3 in Eqs. (jOf), (0) one finds: 



2«o / 

Zi(h) « z -^- M 




1,2. 



(37) 



Thus, Qi ranges from ao to 2ao, or equivalently: zq < Z{(h) < 2zq, on the interval 
of the membrane's curvatures: < 1/r < l/R, or in dimensionless units: < 
h < h . The inverse of Eq. ( |3"TD yields: 



h(zi 



2hQZo 



- 1; 



1,2. 



(3f 



Now, we substitute ( p8| ) into the integral on the right-hand-side of (^6]) and find 
after an integration: 



fi = f l {z i {h),h) 



Pivz Zi(h) 



2hr : 



- 1 



2(1 




(39) 



Substituting Pi = P/2 in (|39D, and then substituting fi from Eq. (|39D into Eq. 
( |10|) instead of fi we find the F(h)/N dependence (in the units of e t ) plotted 
with the dashed curve in Fig. 2. A position of the minimum, defined as the point 
at which dFj dh changes sign, is indicated with the vertical arrow. It is apparent 
from the Fig. 2 that the minimum is rather shallow and happens at the curvature 
radius r comparable with the radius R of the pipette (e.g. of the membrane's 
base). This indicates a substantial deviation of the membrane from the planar 
shape and is in accord with the experiments demonstrating nearly spherical shape 
of the membrane at this pressureaa. Thus, it is not surprising that, unlike in the 
case of a weakly bent thin plateS, now both sides of the membrane's surface are 
stretched. Hence, an important property of this deformed state is the presence in 
both monolayers of nearly equal tensions 7] stretching the membrane and causing 
opening of the protein mechanosensitive channels@i: T\ w T 2 = Pr/A. Let us 
consider another situation, where the symmetry of the tension distribution is 
broken. 



B. Patch-connected case 



Now assume that enough time had elapsed after the beginning of the pressure 
pulse for a transfer of a substantial amount of the lipid molecules from the patch 
to the "sucked in the pipette" part of the membrane (see Fig. 1): t > r. Hence, 
we consider the number of molecules Ni in the patch-connected monolayer as an 
additional variational parameter. In this case the system of equations (p2|), (p3|) 
applies. The first equation written in the dimensionless variables and allowing 
for the general expression (^) yields: 

9 o \ / 1 1 \ ^ 

= 0. (40) 

This equation defines now a zi(h) for the tension-free monolayer with a finite 
curvature ft ^ 0, and thus, substitutes equation (|33]), which is valid for the 
tension-free, but flat case {h = 0). Simultaneously, equation ( |3"TD remains valid 
for the 2-nd monolayer with conserved number of lipids N 2 = No, and provides 
dependence z 2 {h): 






z 2 (h) « z -^ [l 




(41) 



In the Fig. 3 we present our numerical solutions of the Eqs. (^) (dash-dotted 
curve) , (fief ) (solid curves), as well as a dependence given by Eq. (|4]) (dotted 
curve). The insert in Fig. 3 demonstrates that actually z\(h) remains nearly 
equal to zq at all curvatures h. This result is not surprising. Namely, equations 
and (40) differ only by oc h 2 <C 1 and oc h <C 1 terms. This means, that 



the tension-free, relaxed 1-st monolayer, in order to reach its final equilibrium 
state, absorbs as much molecules from the patch as necessary to compensate for 
an increase of its area at a finite curvature. In this way the area per molecule 
remains nearly constant and equal to ~ zq (see insert in Fig. 3). Then, as in 
the flat case, the equilibrium value z\ is again determined by the balance of the 
surface tension and entropic repulsion of the lipid tails. Using dimensionless area 
Zi(h) determined from (|4Ti|), we find the (variable) number of the lipid molecules 
Ni(h), normalized by N , as a function of curvature h, see Fig. 4. This result 
was obtained using the following formula: 



N i =2 Z ° ( h °y I 1 
N Zl (h) \ h) V [h 



h 



2 



(42) 



where equation (|T3|) and definitions ([3(]) where used. 

Now, in order to find the equilibrium value of the curvature h (indicated with 
the arrow in Figs. (3), (4)) at a given pressure difference P, we calculate the total 
free energy F from Eq. (p4"|) normalized by No using the following expression: 



F 
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h{z x {h),h)-^ + f 2 {z 2 {h),h) 
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(43) 



The physical meaning of the factor Ni/Nq in front of fi in combination with 
the last subtracted term oc /o in (|43| ) is simple. Namely, the total number of 
molecules in the system membrane+patch is conserved. Therefore, an increase of 
the number of molecules in the first monolayer, with the free energy per molecule 
fi, occurs at expense of the decrease of the number of molecules in the patch, 
with the energy per molecule Jq. Calculation of the free energy F/Nq, normalized 
by e t , provides its dependence on h expressed with the solid curve in Fig. 2. The 
arrow indicates position of the minimum of the free energy (obtained from the 
calculation of the zero root of the first derivative of the free energy with respect 
to the curvature). 

We have also calculated the dependencies of the dimensionless equilibrium 
curvature h = R/r on the applied pressure difference P for the symmetric strain 
(dashed line): T\ = T 2 , and asymmetric strain (solid line): T\ = 0,T 2 ^ 0, 
equilibrium states, see Fig. 5. As it follows from the figure, the power-law 
dependencies prove to be the same: 1/r ~ const ■ P a , a ~ 1/3, but with the 
constant pre-factor slightly smaller in the symmetric case. 

The influence of the effective surface tension energy f si = 7<2j in Eq. (||), 
characterized by the energy e s (|30D , is demonstrated in Fig. 6. Here, it is obvious, 
that an increase of the surface energy coefficient at a fixed entropic repulsion et, 
causes a greater bending rigidity of the membrane and as a result, a decrease of 
the equilibrium curvature caused by an applied external pressure difference. 

IV. CONCLUSIONS AND EXPERIMENTAL PROPOSALS 

A direct comparison of the two curves in Fig. 2 leads to the main conclusion 
of this work, i.e. the free energy of the relaxed, partially tension-free asymmet- 
ric equilibrium state: T\ = 0, T 2 = Pr/2 is lower than the free energy of the 
non-relaxed symmetric equilibrium state: T\ « T 2 = Pr/4. Thus, the mem- 
brane would change its state until the equilibrium state with the lowest free 
energy would be reached. The number of the lipids in the 1-st monolayer N x in 
the relaxed state is greater than the initial one N by a macroscopical amount 
(~ 15%iV , see Fig. 4). Hence, a process of reaching the relaxed state would 
take time necessary for diffusion of the molecules from the patch into the 1-st 
monolayer through its peripheral edge (see Fig 1). This time, r, as evaluated in 
Section II, Eq. (0), is of order of seconds, and thus is in correspondence with the 
characteristic time of the collapse of the ionic conductance through the protein 
mechanosensitive channels piercing the E-coli's membrane!. Hence, our theory 
confirms the i deJ, that an inter-layer slide in a membrane, soon after the be- 
ginning of the pressure pulse, may cause a relaxation of the lateral stress. This 
stress had initially opened the protein mechanosensitive channels and thus had 



given rise to the ionic conductance. Closing back of the channels within sec- 
onds after the beginning of the pulse would cause simultaneous decrease of the 
ionic conductance as seen experimentallyi. A separate consideration of the two 
equilibrium states discussed above, is justified by the fact that the characteristic 
stress relaxation time: r ~ lsec differs by 5-r6 orders of magnitude from the 
characteristic time At necessary for a propagation of the bending deformations 
over a membrane with the linear dimensions R ~ 1/im: 




V 2 /1 / 3 m-7 \ V 2 

' lg/cm ■ 10 cm x 
10~ 13 erg 



At ~ R 2 — ~ 10 _! W ~ 10~ 5 sec , (44) 



where we used a well known formula for the frequency of the bending waves in a 
thin plateS 

where ks is a bending modulus, p is a mass density, l is a thickness, and q is a 
wave- vector. We also allowed for the role of a bending coefficient ks played by 
the energy et, which follows from direct comparison of the free energy expressions 
©) ® an d the definition of e t given in (|30|). 

Based on our results, shown in Figs. 2, 4, we propose here the following 
possible experimental verifications of the theory. First, it is obvious from the 
Fig. 2 that the equilibrium curvature is greater in the relaxed state than in the 
non-relaxed one, changing correspondingly from 1/r = 0.525-R -1 in the "inter- 
mediate" equilibrium state to 0.675-R -1 in the relaxed state. Hence, we propose 
to measure this change experimentally as an accompanying process of the de- 
crease of the ionic conductance. This change of curvature signifies a simple fact. 
Shortly after the application of the pressure difference P it is balanced by the two 
strained monolayers of a membrane. After enough time has elapsed, a relaxation 
of the strain in the 1-st monolayer occurs, and hence, only the strain in the 2-nd 
monolayer compensates the applied pressure, thus the curvature has to change. 

Second, from Fig. 4 it follows that, as a result of the relaxation of the strain in 
the 1-st monolayer, this monolayer acquires extra lipid molecules from the patch, 
hence the number of molecules in the pipette increases by: (AAq w 15%AT ). As 
this change AiV\ is macroscopic, a corresponding increase of the weight of the 
sucked in membrane during, the decrease of the ionic conductance should be, in 
principle, measurable. 

Finally, we mention some possible future improvements of the theory. These 
could be made by lifting of our simplifying restrictions of the homogeneity of 
the lateral stresses (tensions) across the thickness of each monolayer, as well as 
the approximation of a constant area per molecule over each external surface of 
the membrane. As a consequence, in a more complete scheme of derivations one 
would like not to restrict himself to the spherical symmetry of the membrane's 
shape. Nevertheless, we believe that these improvements can not overturn the 
main physical concept considered in this paper of a diffusion mechanism of a 
stress relaxation in a bilayer membrane with an inter-layer slide. 
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V. FIGURE CAPTIONS 



Fig. 1. Sketch of the experimental setupB. The numbering 1,2 of the mono- 
layers within a bilayer membrane is indicated. 

Fig. 2. Free energy F/(e t N Q ) normalized by e t and number of lipid molecules 
iVo versus dimensionless curvature h. Solid curve corresponds to the asymmetric 
strain state with N2 = Nq < N\. Dashed curve corresponds to the symmetric 
strain case with conserved numbers of lipids in the monolayers: N\ — N2 — Nq. 
Arrows indicate positions of the minima of the free energy. 

Fig. 3. Dimensionless areas per molecule Zi in the i-th monolayer of a bilayer 
membrane as functions of the membrane's curvature h in the asymmetric strain 
state Ti = 0,T 2 = Pr/2 (see text). Arrow indicates equilibrium value of the 
curvature, zq is the equilibrium value under zero external pressure difference. 

Fig. 4. Calculated number of the lipid molecules N± in the first monolayer 
as a function of curvature h in the asymmetric strain state T\ = 0,T 2 = Pr/2. 
Arrow indicates equilibrium value of the curvature at Pv/{2e t ) = 0.0004. 

Fig. 5. Calculated dependencies of the dimensionless equilibrium curvature 
h = R/r on the applied pressure difference P for the symmetric strain (dashed 
line) and asymmetric strain (solid line) equilibrium states in the logarithmic 
scales. 

Fig. 6. Free energy F/(e t N ) normalized by e t and number of lipid molecules 
N Q versus dimensionless curvature h. All notations are as in Fig. 2. The main 
panel corresponds to the ratio t t /t s = 50, which is 10 times greater than for 
the Fig. 2. An insert calculated for the ratio e t /e s = 0.05. Arrows indicate 
equilibrium values of the curvature. 
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